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Abstract 

Objective: Key molecules Involved in microRNA (mlRNA) biogenesis, such as DROSHA, XP05, and DICER, have been 
identified in trophoblast cells, confirming that the miRNA biogenesis pathway is active in human placenta. In addition, 
mlRNAs regulate uterine gene expression associated with inflammatory responses during the peri-implantation period and 
participate in maternal-fetal immune tolerance. The purpose of this study was to demonstrate whether genetic 
polymorphisms in miRNA machinery genes show an association with idiopathic recurrent pregnancy loss (RPL) in Korean 
women. 

Study design: We performed a case-control study with 238 controls and 338 women who had experienced at least two 
consecutive pregnancy losses between 1999 and 2010. Genotypes of miRNA machinery genes, including DICER rs3742330, 
DROSHA rsl0719, RAN GTPase {RAN) rsl4035, and exportin-5 (XP05) rsll077 were analyzed by polymerase chain reaction- 
restriction fragment length polymorphism (PCR-RFLP) assay. The logistic odds ratios (ORs) of RPL were estimated with a 95% 
confidence interval (CI) in multivariate analysis after maternal age adjustment. Gene-gene interactions among the loci of the 
four gene polymorphisms were evaluated using the multifactor dimensionality reduction (MDR) method. 

Resu/ts:The RAN rs14035 CC genotype and DICER rs3742330/D/?OSH/\ rsl0719 GG/TC-nCC, rs3742330//?/A/V rsl4035 GG/CC, 
and DICER rs3742330/XPO5 rsl 1077 GG/AC-nCC combinations were significantly associated with increased RPL risk, whereas 
the RAN rsl4035 CT, DICER rs3742330/R^A/ rsl4035 AA-hAG/CT+TT, DROSHA rs^07^9/RAN rsl4035 TC-nCC/CT+TT, and RAN 
rsl4035/XPO5 rsl 1077 CT-hTT/AA combinations reduced RPL risk. The A-T-T-C and G-C-T-A allele combinations (DICER/ 
DR0SHA/RAN/XP05) were 20 times more frequent in the RPL group than in the control group. 

Conclusion: Our study demonstrates the relationship between RPL development and the polymorphism of the miRNA 
machinery gene RAN and combined genotype of DROSHA/DICER. 
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Introduction 

The incidence of spontaneous pregnancy loss has been 
estimated to be 30% [I]. However, less than 5% of women will 
experience two or more consecutive pregnancy losses [2]. 
Therefore, the loss of two or more clinically recognized 
pregnancies is considered a distinct disorder and defined as 
recurrent pregnancy loss (RPL) for additional evaluation and 
treatment [3]. Despite the various etiologies that result in RPL, 
including cytogenetic abnormalities, antiphospholipid syndrome, 
uterine anomalies, hereditary thrombophilia, autoimmunity, 
sperm quality, and environmental factors, the cause of RPL still 
remains undetermined in the majority of cases. 

MicroRNAs (miRNAs) are ~ 19-25-nucleotide single-strand 
non-coding RNA species that induce post-transcriptional gene 
silencing and mediate translational repression through binding to 



target mRNA, leading to subsequent mRNA degradation. The 
recent elucidation of miRNA function has provided new insight 
into the regulation of gene expression. Key molecules involved in 
miRNA biogenesis, such as DROSHA, XP05, and DICER, have 
been identified in trophoblast cells, confirming that the miRNA 
biogenesis pathway is active in human placenta [4,5]. It has been 
demonstrated that the human placenta produces a large number 
of miRNAs and that those are involved in placental development 
[6]. In addition, miRNAs regulate uterine gene expression 
associated with inflammatory responses during the peri-implanta- 
tion period and participate in maternal-fetal immtme tolerance 
[7,8]. There are numerous reports that demonstrate the associa- 
tion of aberrant miRNA expression with various human diseases 
related to reproductive conditions [9-13]. Several studies have 
reported that genetic polymorphisms are associated with RPL 
development [14-17]. Furthermore, recent studies have found that 
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nucleotide variations witliin pri-miRNA molecules afTect miRNA 
processing and result in altered miRNA expression and that a 
single-nucleotide polymorphism (SNP) in a certain miRNA is 
associated with litter size in pigs [18,19]. 

Although SNPs have been widely implicated in RPL develop- 
ment, such evidence is lacking for miRNA biogenesis pathway 
genes. This study was performed to investigate whether polymor- 
phisms in the miRNA machinery genes DROSHA (rsI0719), 
DICER (rs3742330), RAN (rsl4035), and XP05 (rsll077) are 
associated with the prevalence of RPL in Korean women. Our 
previous study investigated the allele frequencies of nine SNPs in 
those four essential genes in primary ovarian insufficiency and 
control subjects [20]. The results revealed that some polymorphic 
sites, such as XP05 rs3324334 and XP05 rsl 1544382 that had 
been reported to be polymorphic markers, were monomorphic in 
all Korean subjects that we genotyped. Therefore, we examined 
the remaining four polymorphic sites in those four genes. 

Materials and Methods 

Study participants 

The study was reviewed and approved by the institutional 
review board (IRB) of CHA Bundang Medical Center in 1 999, 
and written informed consent was obtained from all participants. 
IRB approved this consent procedure. The study population 
consisted of Korean participants (Asian) recruited from the 
Department of Obstetrics and Gynecology of CHA Bundang 
Medical Center, CHA University (Seongnam-si, South Korea) 
between March 1999 and February 2010. Participants with a 
history of smoking, alcohol abuse, cancer, radiation exposure, 
autoimmune disorder, genetic syndromes, or systemic disease 
affecting ovarian function, such as diabetes mellitus, were excluded 
from the study based on medical history' and physical examination. 
Pregnancy loss was diagnosed with hCG testing, ultrasound, and/ 
or physical examination before 20 weeks of gestational age. RPL 
was defined as two or more consecutive pregnancy losses before 20 
weeks of gestational age. AH participants in the RPL group 
experienced recurrent miscarriage with the same partner. RPL 
patients with previous live births were excluded from this study. 

To diagnose idiopathic RPL, we adopted a set of exclusion 
criteria for the study group. Patients who were diagnosed with 
RPL due to anatomic, chromosomal, hormonal, infectious, 
autoimmune, or thrombotic causes were excluded from the study 
group. To identify anatomic abnormalities in RPL patients, 
sonography, hysterosalpingogram, hysteroscopy, computerized 
tomography, and magnetic resonance imaging were used. 
Karyotyping was performed with standard protocols [21]. 
Hormonal causes included hyperprolactinemia, luteal insufiicien- 
cy, and thyroid disease and were evaluated by measuring levels of 
prolactin, thyroid-stimulating hormone, free T4, follicle-stimulat- 
ing hormone, luteinizing hormone, and progesterone in peripheral 
blood. Lupus anticoagulant and anticardiolipin antibodies were 
examined for autoimmune causes such as lupus and antiphospho- 
Upid syndrome. Thrombotic causes were defined as thrombophilia 
and were evaluated by deficiencies of protein C and protein S and 
by the presence of anti-(32 glycoprotein antibody. The study group 
consisted of 338 women diagnosed with idiopathic RPL [age 
range, 23-43 years; mean age ± standard deviation (SD), 
32.81 ±4.33 years] (Table SI). 

The enrollment criteria for the control group included regular 
menstrual cycles, normal karyotype of 46XX, a history of at least 
one naturally conceived pregnancy, and no history of pregnancy 
loss. The control group was comprised of 238 women [age range, 
22-45 years; mean age ± SD, 33.38±5.79 years]. The average 



gestational age and number of pregnancy losses in the RPL group 
were 7.32±1.85 weeks and 3.04±1.61 losses, respectively. 

Genotyping 

Peripheral blood samples were collected for genotyping. 
Genomic DNA was extracted from collected blood in the presence 
of anticoagulant using a G-DEX blood extraction kit (iNtRON 
Biotechnology, Seongnam-si, Korea). Nucleotide changes were 
determined by polymerase chain reaction-restriction fragment 
length polymorphism (PCR-RFLP) analysis. The PGR primers for 
this study are shown in Table S2. Restriction enzyme digestion 
was carried out using the following enzymes and conditions: jWalll 
{DROSHA rsl 07 19) (New England BioLabs, Ipswich, MA), Banl 
{DICER rs3742330), and EcoRI {XP05 rsl 1077) at 37°C for 
16 hours and BsK (Al/V rsl 4035) at 55°C for 16 hours. Genotypes 
determined by RFLP analysis were confirmed by two independent 
investigators. We confirmed genotypes again by sequencing 10% 
of the samples by random selection. 

Statistical analysis 

The differences in genotype and alkJe combination frequencies 
between idiopathic RPL subjects and controls were compared 
using the multivariate logistic regression and Fisher's exact test, 
respectively. Allele frequencies were calculated to identify devia- 
tions from Hardy-Weinberg equilibrium (HWE) using P= .05 as a 
threshold. Odds ratios (ORs), adjusted odds ratios (AORs), and 
95 '/o confidence intervals (CIs) were used to measure the strength 
of association between genotypes and idiopathic RPL. Two-tailed 
P values <.05 were considered statistically significant. The ORs 
were adjusted by the age of the participants because maternal age 
is a known risk factor for spontaneous abortion and RPL [22]. 

Gene-gene interactions among the SNP loci were analyzed 
using the log-linear model-based multifactor dimensionality 
reduction (LM-MDR) and MDR software version 2.0 (available 
at www.epistasis.org) [23-25]. The allele combination frequencies 
for the selected models by MDR analysis were estimated with the 
HAPSTAT program version 3.0 (www.bios.unc.edu/~lin/ 
hapstat). The false-positive discovery rate (FDR) correction was 
used to adjust multiple comparison tests, and associations with 
FDR-adjusted /"value <.05 were considered to be significant data. 
Statistical analyses were performed with GraphPad Prism 4.0 
(GraphPad Software Inc., San Diego, CA) and StatsDirect 
software version 2.4.4 (StatsDirect Ltd., Altrincham, UK). The 
statistical power was calculated using G*POWER 3.9.1 (Institut 
fur Psychologic, Christian-Albrechts-Universitat Kiel, Kiel, Ger- 
many). 

Results 

The genotypic distribution and allele frequencies of DROSHA 
rsl0719, DICERl rs3742330, A4A'' rsl4035, and XP05 rsl 1077 
are shown in Table 1 . All of the genes that were analyzed showed 
polymorphisms and occurred in HWE in both groups. The DICER 
rs3742330, DROSHA rsl0719, and XP05 rsll077 polymorphisms 
were not associated with the prevalence of RPL. However, the 
Al'V rsl 4035 polymorphism was associated with the prevalence of 
RPL. The frequency of the CT genotype of rsl 4035 was 
significantiy higher in control subjects than in RPL patients. This 
association remained statistically significant after FDR corrc-ction. 
Thus, A4JVrsl4035 CT+TT was associated with the decreased 
prevalence of RPL (AOR = 0.657, 95% CI = 0.469-0.921). 

The genotype frequencies were stratified according to the 
number of consecutive pregnancy losses (Table S3). For DROSHA 
rsl 07 19 and XP05 rsl 1077, there were no statistically significant 
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Table 1. Comparison of genotype frequencies of polymorphisms in miRNA machinery genes between RPL and control subjects. 









RPL patients 










Characteristics 


Controls (n = 238) 


(n = 338) 


AOR (95% CI) 


P' 


Statistical power (%) 


DICER rs3742330 


AA 


75 (31.5) 


119 (35.2) 


1 .000 (reference) 








AG 


123 (51.7) 


1 52 (45.0) 


0.741 (0.510-1.077) 


0.117 


0.234 


42.5 


GG 


35 (14.7) 


67 (19.8) 


1.196 (0.724-1.975) 


0.484 


0.776 


12.5 


Dominant (AA vs. AG+GG) 






0.841 (0.590-1.198) 


0.336 


0.672 


15.4 


Recessive (AA+AG vs. GG) 






1.428 (0.913-2.235) 


0.119 


0.476 


48.4 


HWE P 


0.18 


0.147 










DROSHA rs10719 


TT 


110 (46.2) 


161 (47.6) 


1 .000 (reference) 








TC 


108 (45.4) 


1 50 (44.4) 


0.943 (0.667-1.335) 


0.742 


0.977 


6.4 


CC 


20 (8.4) 


27 (8.0) 


0.919 (0.487-1.733) 


0.794 


0.794 


5.6 


Dominant (TT vs. TC+CC) 






0.933 (0.669-1.303) 


0.685 


0.859 


7.0 


Recessive (TT+TC vs. CC) 






0.907 (0.494-1.665) 


0.752 


0.954 


6.4 


HWE P 


0.363 


0.329 










RANts^'^o^5 


CC 


123 (51.7) 


210 (62.1) 


1 .000 (reference) 








CT 


104 (43.7) 


1 1 3 (33.4) 


0.640 (0.452-0.906) 


0.012 


0.048 


72.6 


TT 


11 (4.6) 


15 (4.4) 


0.796 (0.354-1.791) 


0.582 


0.776 


72.5 


Dominant (CC vs. CT+TT) 






0.657 (0.469-0.921) 


0.015 


0.060 


70.0 


Recessive (CC+CT vs. TT) 






0.977 (0.439-2.1 72) 


0.954 


0.954 


4.6 


HWE P 


0.059 


0.967 










xposts^^o^^ 


AA 


197 (82.8) 


279 (82.5) 


1 .000 (reference) 








AC 


39 (16.4) 


53 (15.7) 


0.993 (0.632-1.561) 


0.977 


0.977 


50.5 


CC 


2 (0.8) 


6 (1.8) 


1.896 (0.357-10.061) 


0.453 


0.776 


12.6 


Dominant (AA vs. AC+CC) 






1.041 (0.670-1.617) 


0.859 


0.859 


53.4 


Recessive (AA+AC vs. CC) 






1.911 (0.360-10.140) 


0.447 


0.894 


2.7 


HWE P 


0.964 


0.072 











Note: For AOR, OR was adjusted by age of participants. RPL = recurrent pregnancy loss; AOR = adjusted odds ratio; CI = confidence interval; HWE = Hardy-Weinberg 

equilibrium; 

^Fisher's exact test; 

''FDR-adjusted P value. 

doi:1 0.1 371 /journa!.pone.0095803.t001 



differences between the two groups. The DICER rs3742330 
recessive model (GG compared with AA+AG) was associated with 
subjects wlio had history of three or more and four or more 
pregnancy losses (AOR= 1.718, 95% 01= 1.032-2.861 and 
AOR = 2.170, 95% 01= 1.181-3.986, respectively). For RAM 
rs 14035, the OT genotype was less frequently observed in the 
subgroups of RPL participants with two or more and four or more 
consecutive pregnancy losses (AOR = 0.624, 95% 01 = 0.412- 
0.947 and AOR =0.502, 95% 01 = 0.287-0.879, respectively). 

We performed combination analyses for miRNA machinery 
gene polymorphisms that we tested for the study (Table 2). These 
analyses revealed that the DICER rs3742330 GG/DROSHA 
rsl0719 TO+00 genotype was significantly more frequent in 
RPL patients than in control subjects (AOR = 1.990, 95% 
01= 1.048-3.778). The frequency of the DICER rs3742330 AA+ 
AG/AWrs 14035 OT+TT genotype was significantly higher in 
the control group (AOR = 0.688, 95% 01 = 0.475-0.997). In 
addition, DROSHA rsl0719 TO+OO/AIA' rsl4035 OT+TT and 
iLLA''rsl4035 0T+TT/ZP05 rsll077 AA genotypes were more 



frequently observed in control participants (AOR = 0.621, 95% 
01 = 0.389-0.991, AOR = 0.631, 95% 01 = 0.435-0.914, respec- 
tively). Oombination analysis according to the number of 
pregnancy losses was also performed (Table 3). Although these 
associations with RPL were not significant after FDR correction in 
combination analysis (Table 2), those associations became more 
evident when we performed combination analysis according to the 
number of pregnancy losses even after FDR correction (Table 3). 
In this combination analysis, there was a significant association 
between RPL development and the DICER/DROSHA GG/TO+ 
00, DICER/RAM GG/00, and DICER/XP05 GG/AO+00 
genotypes. The RAM OT+TT genotype showed a significant 
association with RPL when it was combined with the XP05 AA 
genotype. These significances were maintained when we per- 
formed multiple comparison correction using FDR, for all except 
the RAN/XP05 OT+TT/AA genotype. Thus, both individual and 
combined genotyping analysis suggested that both the DICER 
rs3742330 GG and RAN rsl4035 00 genotypes are associated 
with RPL prevalence. 
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Table 2. Combination analysis of polymorphisms of miRNA machinery genes in RPL patients and control participants. 









RPL patients 








Statistical power 


Genotypes 


Controls (n = 


238)(n = 338) 


AOR (95% CI) 




(%) 


DICER/DROSHA 


AA+AG/TT 


90 (37.8) 


141 (41.7) 


1 .000 (reference) 








AA+AG/TC+CC 


113 (47.5) 


130 (38.5) 


0.729 (0.505-1.052) 


0.091 


0.273 


45.7 


GG/TT 


20 (8.4) 


20 (5.9) 


0.633 (0.322-1.245) 


0.185 


0.370 


25.5 


GG/TC+CC 


15 (6.3) 


47 (13.9) 


1.990 (1.048-3.778) 


0.036 


0.138 


64.0 


DICER/RAN 


AA+AG/CC 


107 (45.0) 


168 (49.7) 


1 .000 (reference) 








AA+AG/CT+TT 


96 (40.3) 


103 (30.5) 


0.688 (0.475-0.997) 


0.048 


0.273 


56.5 


GG/CC 


1 6 (6.7) 


42 (12.4) 


1.688 (0.903-3.156) 


0.101 


0.303 


42.0 


GG/CT+TT 


1 9 (8.0) 


25 (7.4) 


0.826 (0.433-1.576) 


0.562 


0.674 


8.4 


DICER/XP05 


AA+AG/AA 


1 67 (70.2) 


226 (66.9) 


1 .000 (reference) 








AA+AG/AC+CC 


36 (15.1) 


45 (13.3) 


0.963 (0.592-1.566) 


0.878 


0.905 


5.3 


GG/AA 


30 (12.6) 


53 (15.7) 


1.308 (0.801-2.137) 


0.283 


0.425 


20.4 


GG/AC+CC 


5 (2.1) 


14 (4.1) 


2.055 (0.725-5.822) 


0.175 


0.350 


32.9 




TT/CC 


60 (25.2) 


102 (30.2) 


1 .000 (reference) 








TT/CT+TT 


50 (21.0) 


59 (17.5) 


0.694 (0.424-1.137) 


0.147 


0.294 


37.6 


TC+CC/CC 


63 (26.5) 


108 (32.0) 


1.024 (0.654-1.604) 


0.917 


0.917 


5.2 


TC+CC/CT+TT 


65 (27.3) 


69 (20.4) 


0.621 (0.389-0.991) 


0.046 


0.138 


65.0 




TT/AA 


90 (37.8) 


133 (39.3) 


1 .000 (reference) 








TT/AC+CC 


20 (8.4) 


28 (8.3) 


0.945 (0.502-1.781) 


0.861 


0.905 


5.3 


TC+CC/AA 


107 (45.0) 


146 (43.2) 


0.914 (0.633-1.319) 


0.631 


0.757 


73.3 


TC+CC/AC+CC 


21 (8.8) 


31 (9.2) 


1.016 (0.548-1.886) 


0.959 


0.959 


5.0 




CC/AA 


103 (43.3) 


1 77 (52.4) 


1 .000 (reference) 








CC/AC+CC 


20 (8.4) 


33 (9.8) 


0.964 (0.525-1.768) 


0.905 


0.905 


5.1 


CT+TT/AA 


94 (39.5) 


102 (30.2) 


0.631 (0.435-0.914) 


0.015 


0.090 


73.3 


T+TT/AC+CC 


21 (8.8) 


26 (7.7) 


0.708 (0.376-1.334) 


0.285 


0.428 


17.7 



Note: For AOR, OR was adjusted by age of participants. RPL = recurrent pregnancy loss; AOR = adjusted odds ratio; CI = confidence interval. 

^Fisher's exact test; 

''FDR-adjusted P value. 

doi:1 0.1 371 /journal.pone.0095803.t002 



To evaluate whether gene-gene interaction models exert 
synergistic effects on RPL risk, we conducted MDR-based allele 
combinations of four polymorphisms of miRNA machinery genes 
(Table 4 and Table S4). We selected two interaction models 
{DICER/ DROSHA/ RAN/ XP05 and DICER/DROSHA) because 
their CV values were 10/10. Among the models of four 
polymorphic loci in DICER/ DROSHA/ RAN/ XP05, four aUele 
combinations (A-T-C-A, A-T-T-C, G-T-C-C, and G-C-T-A) 
showed synergistic effects on increased RPL risk (OR =1.510, 
95% CI= 1.161-1.965; OR=20.858, 95% 01= 1.240-350.756; 
OR= 12.117, 95% CI = 0.697-2 10.595; and OR = 29.759, 95% 
CI= 1.794-493.586, respectively), while A-T-C-C, A-C-T-A, and 
G-T-T-A were associated with reduced RPL risk (OR = 0.366, 
95% CI = 0.189-0.709; OR=0.466, 95% CI = 0.281-0.773; and 
OR =0.533, 95% CI = 0.335-0.850, respectively). However, it 
was difficult to make a conclusion with four allele combinations 
{DICER/ DROSHA/ RAN/ XPOj) of A-T-T-C, G-T-C-C, and G-C- 
T-A due to the wide range of confidence intervals. Among the 



models of two loci, the A-T and G-C allele combinations in the 
DICER/DROSHA model were associated with increased preva- 
lence of RPL (OR= 1.292, 95% CI= 1.017-1.643 and 
OR= 1.694, 95% CI = 1.182-2.426, respectively), whereas A-C 
was associated with reduced RPL risk (OR= 0.615, 95% 
CI = 0.451-0.839). Table S5 represents the statistical power of 
the significant genotype. In addition, we conducted power analysis 
to detect the associations in the case-control study at the 5% 
significance level. The power varies for an OR by the proportion 
of exposure (i.e., the target genotype) in the control. A total sample 
size of 576 (338 RPL cases vs. 238 controls) could have reasonable 
power to detect an OR of 1 .7 when the proportions of exposure in 
the control were 20% (power = 76.3%), 30% (power = 84.9%), and 
40% (power = 87.7%). Additionally, this sample size has reason- 
able power to detect an OR of 0.55 for the proportions of 20% 
(power =73.1%), 30% (power = 85.6%), and 40% (pow- 
er =91.3%) in the control. Furthermore, this sample size might 
have sufficient power to detect an OR of 1 . 7 or greater and an OR 
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Table 4. Allele connbination analysis of mlRNA machinery genes In RPL patients and control participants using the multifactor 
dimensionality reduction method. 





Allele combination 


Controls (n = 238) 


RPL patients 
(n = 338) 


OR (95% CI) 






Statistical power (%) 


DICER/DR0SHA/RAN/XP05 


A-T-C-A 


0.2455 


0.33 


1.510 (1.161-1.965) 


0.002 


0.010 


71.4 


A-T-C-C 


0.0548 


0.0211 


0.366 (0.189-0.709) 


0.003 


0.010 


70.0 


A-T-T-A 


0.0726 


0.0636 


0.856 (0.539-1.359) 


0.552 


0.589 


11.5 


A-T-T-C 


0.0001 


0.021 


20.858 (1.240-350.756) 


0.001 


0.008 


8.6 


A-C-C-A 


0.1175 


0.0874 


0.717 (0.487-1.055) 


0.110 


0.176 


31.7 


A-C-C-C 


0.0057 


0.0143 


2.367 (0.648-8.651) 


0.259 


0.345 


26.8 


A-C-T-A 


0.0812 


0.0395 


0.466 (0.281-0.773) 


0.003 


0.010 


66.6 


A-C-T-C 


0.0067 


0 


0.100 (0.005-1.942) 


0.070 


0.140 


16.6 


G-T-C-A 


0.2176 


0.195 


0.868 (0.650-1.159) 


0.336 


0.384 


16.4 


G-T-C-C 


0 


0.0114 


12.117 (0.697-210.595) 


0.024 


0.055 




G-T-T-A 


0.0893 


0.0504 


0.533 (0.335-0.850) 


0.008 


0.021 


59.4 


G-T-T-C 


0.0092 


0.0057 


0.702 (0.175-2.823) 


0.724 


0.724 


12.3 


G-C-C-A 


0.086 


0.1077 


1.284 (0.859-1.920) 


0.231 


0.336 


21.5 


G-C-C-C 


0.0083 


0.0215 


2.678 (0.883-8.121) 


0.098 


0.174 


37.0 


G-C-T-A 


0 


0.0302 


29.759 (1.794-493.586) 


<.0001 


0.002 




G-C-T-C 


0.0056 


0.001 1 


0.234 (0.024-2.254) 


0.313 


0.384 


9.1 


DICER/DROSHA 


A-T 


0.3761 


0.4373 


1.292 (1.017-1.643) 


0.039 


0.052 


43.8 


A-C 


0.208 


0.1397 


0.615 (0.451-0.839) 


0.002 


0.008 


69.3 


G-T 


0.313 


0.261 


0.773 (0.596-1.001) 


0.054 


0.054 


38.9 


G-C 


0.1029 


0.1621 


1.694 (1.182-2.426) 


0.004 


0.008 


67.9 


Note: ORs and 95% CIs of each allele combinations were calculated with reference to frequencies of all others 
RPL = recurrent pregnancy loss; AOR = adjusted odds ratio; CI = confidence interval. 


using Fisher's exact test. P value by Fisher's exact test. 



^Fisher's exact test; 
"FDR-adjusted P Value. 
doi:1 0.1 371 /journal.pone.0095803.t004 



of 0.55 or less. However, the power will be smaller when it is 
compared with multiple genotypes. 

Discussion 

We evaluated four SNPs in DICER, DROSHA, XP05, and RAN, 
which are essential for miRNA biogenesis [26,27]. Our results 
comparing differences in genotype frequencies in miRNA 
machinery genes revealed that DICER rs3742330 GG genotype 
had significant associations with RPL in Korean women when it 
was combined with the DROSHA rsl0719 TC+CC genotype. In 
addition, the ORs of RPL prevalence were higher because patients 
who possessed the DICER rs3742330 GG/DROSHA rsl0719 TC+ 
CC genotype experienced more pregnancy loss. Defective DICER 
function due to genetic polymorphism might be intensified by 
DROSHA polymorphisms. These two molecules are considered 
major regulators of miRNA biogenesis. However, there is no 
evidence that these molecules interact with each other. It remains 
unclear how genetic polymorphisms in these genes influence each 
other. To elucidate the mechanisms by which different alleles of 
DICER or DROSHA affect RPL development, additional function- 
al studies should be performed. Although the exact mechanism 
could not be elucidated in the present study, our finding might be 
used in the future as a marker to predict RPL development in 
women with a spontaneous abortion. 



Dicer is essential in miRNA biogenesis, and its function has 
been reported to be involved reproduction. Dicer expression is 
detected in female reproductive organs, including the ovary, 
oviduct, and uterus [28-30]. Female mice with general hypomor- 
phic mutation of Dicer were infertile because of luteal deficiency 
[30]. The conditional knockout oi Dicer within ovarian granulosa 
cells affects ovulation rate by influencing the number of 
preovulatory foUicles that achieve proper development [28,29]. 
On the other hand. Dicer is reported to be involved in regulating 
placental growth by influencing pro-mitogenic signaling molecules 
through miRNA action [31]. 

Although there is limited evidence that polymorphisms of 
DICER alter biologic functions depending on alleles, a number of 
studies have reported that DICER polymorphisms, including 
rs3742330, affect disease development and patient survival in 
various cancers [31-33]. In addition, the functional analysis of the 
DICER rsl057035 SNP, which resides in the 3'UTR oi DICER, 
has revealed that targeting of has-miR-574-3p to DICER was 
differentially affected by the type oi DICER rs 105 7035 allele and 
that the DICER rsl057035 variant C aUele led to reduced 
expression of Dicer [34] . In the DICER gene, the polymorphic site 
DICER rs3742330A>G, found in the 3' UTR of DICER, is 
considered to be important for mRNA transcript stability. In terms 
of human reproduction, DICER rs 12323635, located in promoter 
region of DICER, has been reported to decrease the risk of 
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oligozoospermia [35]. Therefore, Dicer is likely to be a crucial 
molecule for pregnancy. 

There are several potential mechanisms that explain why DICER 
polymorphisms are associated with miscarriage. DICER polymor- 
phism may affect decidualization of the endometrium and induce 
implantation defects in embryos, resulting in recurrent miscarriage. 
Decidualization of endometrial stroma is one of the crucial steps for 
trophoblast invasion and placenta development [36] . A recent study 
using human embryonic stem cells demonstrated that DICER 
expression increased during in vitro decidualization. DICER 
inhibition affected the expression of decidual markers and 
transcription factors — such as PRE and HOXAIO — which play a 
role during decidualization [37]. Second, DICER polymorphisms 
may have an influence on placental development by affecting 
miRNA expression involving cytotrophoblast proliferation. Cyto- 
trophoblasts are trophoblastic stem cells that play multiple roles [6] . 
Recent research has demonstrated enhanced proliferation of 
cytotrophoblasts without exogenous growth factors following Dicer 
knock-down in first-trimester placenta. Although the miRNAs 
affecting cell proliferation were not identified, the authors suggested 
that Dicer-dependent miRNA affected the proliferation of cytotro- 
phoblasts [5] . Dicer polymorphisms induce aberrant expression of 
miRNAs regulating cytotrophoblast proliferation, possibly resulting 
in abnormal placental development and abortion. Finally, the 
miRNA expression affected by the DICER polymorphism may 
induce chromosome mis-segregation, leading to aneuploidy, which 
is the most common cause of miscarriage. The centrosome is known 
to be a source of aneuploidy. The extra centrosome forms a 
multipolar spindle during mitosis in a mother cell, resulting in 
aneuploid daughter cells. This pathogenesis has been shown in 
cancer development [38]. miR-210 overexpression is reported to 
induce centrosome amplification in renal cancer [39]. Therefore, 
aberrant miRNA expression may cause chromosomal aneuploidy in 
cancer. A similar mechanism exists by which abnormal miRNA 
expression may cause aneuploid embryos and recurrent miscarriage 
during early development. A significant association between the 
Al'Vrsl4035 genotype and RPL prevalence was observed when it 
was combined with the XPOJrsl 1077 AA genotype. Ran, which is a 
member of the Ras superfamily of GTPases, is essential for the 
translocation of prc'-miRNAs from the nucleus to the cytoplasm 
through the nuclear pore complex [40]. Xpo5 mediates nuclear 
export of pre-miRNA in a RAN GTP-dependent manner by 
binding to pre-miRNA and RAN GTPase in the nucleus [40]. 
Although several reports have shown the association of these 
molecules with disease development in various cancers, no report 
has elucidated the role of these mc)lecules in human reproduction 
[41,42]. 

Our previous study revealed that the XP05 rs2257082 allele was 
more frequently observed in premature ovarian insufficiency (POT) 
patients than in control subjects [20]. However, we could not 
identify any association between RAJ^ polymorphism and POI 
development. Additionally, iM-Vrs 1403,5, but not XP05 rsll077, 
exerted an effect on RPE development. Regarding both molecules, 
few reports have evaluated the clinical significance of SNPs not 
only in cancer but also in pregnancy-related disease. In addition, 
no study has elucidated the functions of those SNPs in miRNA 
biogenesis. Therefore, additional functional studies are required to 
elucidate how RAN rs 14035 affects disease promotion and 
prevention depending on its genotype. 

We obtained contrasting results regarding RPL development in 
tiie combination analysis for the DROSHA rsl0719 TC-I-CC 
genotype. When it was combined with DICER rs3724330 GG 
genotype, RPE development was more frequentiy observed in 
RPE group. When die DROSHA rsl0719 TC-hCC genotype was 



combined with the RAN rs 14035 CT-I-TT genotype, RPE 
development was rare in the RPE group. One potential 
explanation for this observation is that idiopathic RPE has various 

unknown etiologies. Many reports in the literature have demon- 
strated that several causes are associated with recurrent miscar- 
riage, such as angiogenesis, immunologic factors, and miRNAs 
[15,18,43]. This may explain why the same polymorphism affects 
disease development differently. The result of our combination 
analysis of polymorphisms according to the number of pregnancy 
losses supports this explanation (Table 3). The AORs increased as 
the number of previous pregnancy losses increased in patients who 
have the DROSHA rsl0719 TC-I-CC genotype combined witii tiie 
DICER rs3742330 GG genotype (1.596 vs. 2.460 vs. 4.084). 
Conversely, no such tendency was noted with the DROSHA TC+ 
CC/AiA''CTH-TT genotype. Instead, the effect of polymorphisms 
was attcnuat(;d as the number of miscarriages increased. This 
suggests that genetic polymorphisms with respect to miRNA 
machinery genes were associated with disease development not 
through a common pathway but through various pathogeneses. 

There are several potential limitations in this study. First, this 
case-control study only shows an association between SNPs in 
miRNA machinery genes and RPE development, not a causal 
relationship. Therefore, it is difficult to conclude that miRNA 
machinery genes affect recurrent miscarriage. Second, the manner 
in which polymorphisms in miRNA biogenesis pathway genes 
affect RPE development is still unclear, and a functional study of 
those SNPs to elucidate the pathogenesis related to RPE was not 
conducted. Third, the lack of information regarding placental 
pathology, immunologic profiles, and miRNA expression could 
have contributed to the investigation of potential roles of miRNA 
machinery genes during the peri-implantation and early pregnan- 
cy periods. Eastiy, the sample size of the control group is relatively 
small, compared with that of the study group. This is because we 
included all eligible cases during the study period to avoid selection 
bias. This difference in sample size wiU increase type 2 errors in 
the genetic association study. To reduce the bias, we compared the 
AW rs 14035 and XP05 rsll077 genotype frequencies in the 
control group with those from other studies that were performed in 
Korea. The genotype frequencies of RAM rs 14035 and XP05 
rsl 1077 in the control group of our study showed similar patterns 
to the results in the control groups of previous reports regarding 
lung cancer and breast cancer in Korea [26,27]. 

Recent groundbreaking discoveries that small, non-coding 
RNAs can regulate gene expression transcriptionally or post- 
transcriptionally and that these RNAs are involved in placental 
development through trophoblast cUfferentiation, proliferation, 
and angiogenesis provide us with insight into the pathogenesis of 
idiopathic RPE [4,7,44]. To the best of our knowledge, this is the 
first report to investigate the association between polymorphisms 
in miRNA biogenesis machinery genes and the ])rc\ alcnce of RPE 
in a Korean population and suggests that genetic polymorphisms 
in miRNA biogenesis machinery genes show a relationship with 
the risk of RPE. The present study may contribute to assessments 
of individual risk of RPE. However, further epidemiologic studies 
with larger subject numbers should be performed to confirm and 
expand our results. In addition, the results of our study warrant 
further functional studies to elucidate the mechanisms by which 
polymorphisms of miRNA machinery genes affect RPE develop- 
ment. 
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